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Abstract: We perform a detailed analysis of the SO(IO) SUSY GUT model with D3 family 
symmetry of Dermisek and Raby (DR). The model is specified in terms of 24 parameters 
and predicts, as a function of them, the whole MSSM set of parameters at low energy scales. 
Concerning the SM subset of such parameters, the model is able to give a satisfactory 
description of the quark and lepton masses, of the PMNS matrix and of the CKM matrix. 
We perform a global fit to the model, including flavour changing neutral current (FCNC) 
processes Bs — > fJ'~^fJ'~, B Xg^, B Xsi~^£~ and the B^^s — Bd,s mass differences AM^^c, 
as well as the flavour changing (FC) process B^ — s- r^z^. These observables provide at 
present the most sensitive probe of the SUSY mass spectrum and couplings predicted by the 
model. Our analysis demonstrates that the simultaneous description of the FC observables 
in question represents a serious challenge for the DR model, unless the masses of the 
scalars are moved to regions which are problematic from the point of view of naturalness 
and probably beyond the reach of the LHC. We emphasize that this problem could be 
a general feature of SUSY GUT models with third generation Yukawa unification and 
weak-scale minimal flavour violation. 
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1. Introduction 

A well-known problem of supersymmetry (SUSY) at the electroweak (EW) scale is its 
proliferation of parameters, arising if one keeps the most general allowed terms in the soft 
sector. In absence of the latter, SUSY is not phenomenologically viable, while in presence 
of a most general soft sector, SUSY largely loses its predictivity, due to the bulkiness of 
the parameter space. 
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Since SUSY goes often together with the idea of Grand Unification, due to the tanta- 
lizing observation of gauge coupling unification after MSSM running, a sensible way-out to 
the above problem is to take a 'top-down' approach. In this case one starts with the theory 
at the GUT scale - with a far simpler parameter space than the MSSM one - and runs 
all the parameters downwards through Renormalization Group Equations (RGEs). The 
latter then dynamically generate all the mass splittings for the soft terms at the EW scale, 
which, in the low-energy MSSM without GUT, are treated instead as free parameters. 

As an argument in favour of a top-down approach to the MSSM, it can be noted that 
predictive SUSY GUT models need typically two ingredients at the GUT scale. The first 
one is the choice of a specific GUT gauge symmetry. Second, one has to fix the high-scale 
soft sector with minimal assumptions on the parametric dependence, which seems justified 
due to a higher amount of symmetry of the theory at this high scale. A further restriction 
on the number of parameters and hence more predictability of the model can be obtained 
by introducing additional family symmetries. Once these different parts of the model are 
specified, the theoretical and computational tools available today allow for a very controlled 
theoretical error on the model parameters at low-energy, in spite of the 'long running' from 
the GUT scale and the presence of mass thresholds. 

Concerning SUSY GUT models present on the market, while it is easy to construct 
models reproducing the gauge sector of the MSSM at low-energy, it is far more challenging 
to find models also correctly describing 'flavour patterns' such as quark and lepton masses 
and the CKM and PMNS matrices. 

One such notable model has been proposed by Dermisek and Raby in It is an 
SO(IO) SUSY GUT, augmented with a simple family symmetry, conserving R-parity at 
low-energy. As shown in , this model is able to successfully fit all the parameters of the 
SM. In particular, by using some observables basically unaffected by SUSY contributions, 
it can reproduce the CKM matrix entries. Finally, it also describes the known parameters 
in the neutrino sector. In Ref. [^, the same model was also extensively studied in the 
lepton flavour sector, thereby providing a number of signatures of the model, which should 
at least in part be tested by forthcoming experiments. 

In view of the above mentioned remarkable performance of the DR model in describing 
low-energy observables, it is interesting to have a closer look at its SUSY spectrum, with 
the aim of testing, e.g., the predicted mass hierarchies. Since SUSY particles have not yet 
been observed, such task can only be accomplished by analyzing implied loop effects in 
flavour changing neutral current (FCNC) processes. In fact, the specific mass patterns of 
gauginos, up- and down-squarks and Higgs multiplets predicted by the model, do affect 
measured quark FCNC observables in a peculiar way, and if such observables are well 
controlled, the pattern of implied corrections is testable. Fortunately, we have today a 
whole host of such observables, which have the virtue of being at the same time precisely 
measured and accurately calculated within the MSSM. 

The aim of the present paper is then an in-depth test of the model in the light of all 
the most accurate information presently available on quark FCNCs. The flavour sector 
is often overlooked in first analyses of new physics (NP) models, due to its vastness and 
the necessity of sometimes involved calculations. In the context of the present model, we 
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show that the flavour sector has nonetheless enough sensitivity to the details of the SUSY 
spectrum, to represent a discriminating test for the model itself. 

Our paper is organized as follows. In Section ^ we give a brief overview of the DR 
model, focusing on those ingredients that are most relevant for our purposes. In Section ^ 
we then outline the procedure to connect the GUT scale model with low-energy observables. 
Such procedure is well-known, but often obscured by the assumptions made on the running 
and the mass thresholds. Section ^ goes then in more detail on the determination of masses 
and couplings, through consideration of low-energy threshold effects. In Section ^ we then 
present the collection of FCNC observables we use for the purpose of our paper. The 
emphasis here is on presenting simplified expressions, in order to provide an intuitive picture 
of the main effects, with refined formulae only used in the actual numerical analysis. In 
the light of such intuitive expressions. Section ^ then discusses the general pattern featured 
by the corresponding FCNC observables within the DR model. An extensive analysis of 
all such features upon variation of the DR model parameters is then presented in Section 
1^. Finally, Section ^ is devoted to our conclusions. 

2. The model 

The DR model |l|, ^ is a supersymmetric SO(IO) Grand Unified Theory with an additional 
D3 X [U{1) X Z2 X Z3] family symmetry. 

The above symmetry group fixes the following structure for the superpotential 

W = Wf + , 

with 

Wf = 163l0163 + 16alOXa 

+ Xa{Mx Xa + 45 ^ I63 + 45 t 16a + A 16„) , (2.1) 

MM 

= T6(A2 Na 16a + A3 iV3 I63) + ^iSa Na Na + ^3 N3 N3) . (2.2) 

The first two families of quarks and leptons are contained in the superfleld IGq, a = 
1,2, which transforms under SO(10)xD3 as (16,2^), whereas the third family in I63 
transforms as (16, 1^). The two MSSM Higgs doublets and Hfi are contained in a 
10. As can be seen from the first term on the right-hand side of (p.l|), Yukawa unification 
M = \ = = ^u-r at Mg is obtained only for the third generation, which is directly 
coupled to the Higgs 10 representation. This immediately implies large tan /3 ~ 50 at low 
energies. 

The effective Yukawa couplings of the first and second generation fermions are instead 
generated via the Froggatt-Nielsen mechanism Q as follows. Additional fields are intro- 
duced, i.e. the 45 which is an adjoint of SO(IO), the SO(IO) singlet fiavon fields (j)"" , (j)"- , A 
and the Froggatt-Nielsen states Xa,Xa- The latter transform as a (16,2^) and a (16,2^), 
respectively, and receive masses of 0{Mg) as acquires an SO(IO) breaking VEV. Once 



-3- 



they are integrated out, they give rise to effective mass operators which, together with the 
VEVs of the ffavon fields, create the Yukawa couplings for the first two generations. This 
mechanism breaks systematically the full flavour symmetry and produces the right mass 
hierarchies among the fermions. 

The obtained Yukawa matrices for up-quarks, down-quarks, charged leptons and neu- 
trinos are 




Y, 



Y. 




(2.3) 



From eqs. (^]^) one can see that the flavour hierarchies in the Yukawa couplings are encoded 
in terms of the four complex parameters /?, a, £, ^ and the additional real ones e,e',A. 

In order to avoid neutrino masses of the order of the other fermion masses, one invokes 
the type-I see-saw mechanism |5|, ^, 0, |8|. In particular, three SO(IO) singlet Majorana fer- 
mion fields Na, {a = 1^2) are introduced via the contribution of ^ {Sa Na Na + S3 N^) 
to the superpotential. 

The mass term ^ N is produced when the flavon fields acquire VEVs (Sa) = 

and {S3) = Matj. Together with a 16 Higgs one is allowed to introduce the interaction 
terms 16 (A2 Na 16^ -|- A3 N3 I63), which in turn generate a mixing matrix V between the 
right-handed neutrinos and the additional singlets (y V N), when the 16 acquires an SO(IO) 
breaking VEV (16)f = ^^le- The resulting effective right-handed neutrino mass terms read 



Wn 



VV N + -N Mn N , 



(2.4) 



/ A2 



V 



V16 



A2 
A3 



Mn = diag(M^, , Mn^ , Mn, 



Diagonalization leads to the effective right-handed neutrino Majorana mass 
Mr = -V = -diag(M^j, , Mr, ,Mr,) . 



(2.5) 



(2.6) 



By integrating out the EW singlets V and N, which both receive GUT scale masses, one 
ends up with the light neutrino mass matrix at the EW scale given by the usual see-saw 
formula 



M = 1711, Mo ^ ml . 



(2.7) 



3. Basic procedure 

In this section, we describe the procedure which, from the specification of the model at the 
GUT scale, leads to the MSSM mass spectrum and observables at the EW scale. 
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3.1 Step 1: Parameters of the model 
Gauge coupling sector 

We choose as three parameters the unification scale Mq, the gauge coupling defined 
through 

ac = aiiMc) = a2{MG) , (3.1) 
and the threshold correction es defined through 

aslMc) = 0^(1 + 63) . (3.2) 
The threshold correction 63 helps to obtain the right value of 03 (M^). 
SUSY sector 

We have the following set of soft SUSY breaking parameters: a universal sfermion mass 
mi6, a universal gaugino mass Mi a universal trilinear coupling parameter Aq, and the 
Higgs mass parameters mu^ and mriH^- It is well-known that at large tan /?, EWSB is easier 
to achieve by allowing the soft SUSY breaking Higgs mass parameters to be split. This 
also has the consequence that the absolute value of ^ is not fixed by EWSB, as is the case 
in the CMSSM^, but is instead a free parameter. 

Yukawa matrices 

The Yukawa matrices for up- and down- type quarks, charged leptons and neutrinos are 
parameterized as given in (plD , with e,e',A being real and p, a, e, ^ complex parameters. 



Right-handed neutrinos 

The diagonal right-handed neutrino mass matrix is chosen to be real, which amounts to 
three parameters M^^, with i = 1,2,3. 

Weak-scale parameters 

In addition to the above GUT-scale parameters, the SUSY parameters p and tan /? have 
also to be specified as an input at the weak scale. 

The total number of input parameters, listed in this step and summarized in Table |l|, 
is therefore 24. Once these parameters are fixed, the MSSM couplings as well as its whole 
mass spectrum (which includes the SM part) can be predicted at energies below Mq, in 
particular at Mz or any lower scale relevant for FCNC processes. Of course a subset of 
the observables amenable to prediction is used to fix the model parameters. 

The procedure to fix the above parameters will be addressed in detail in Section |^. 
Such procedure uses low-energy observables, and for this reason we need to evolve the 
fundamental parameters of the theory to Mz or below, using their RGEs. The evolution 
to low energy is performed according to a procedure we will now describe. 



^Very interesting combined analyses in the framework of the CMSSM have been performed in j^, |l^, p^f . 
In jl^], also non-universal Higgs masses were considered, but only with vanishing Aq. 
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Sector 


# 


Parameters 


gauge 


3 


"G, Mg, €3, 


SUSY (GUT scale) 


5 


mid, Ml/2, Ao, mu^, mH^, 


textures 


11 


e, e', X, p, a, e, i. 


neutrino 


3 


Mr,, Mr,, Mr,, 


SUSY (EW scale) 


2 


tan P, fi 



Table 1: Parameters in the DR model. 



3.2 Step 2: RG evolution 
Mgut > scale > Mr, 



The RGEs of the MSSM with right-handed neutrinos |]l^, 14, 15| are used to run the gauge 
couplings, Yukawa couplings, all soft SUSY breaking parameters and the right-handed 
neutrino mass matrix Mr down to the mass of the heaviest right-handed neutrino. At 
this threshold, one must rediagonalize Mr and integrate out one neutrino following ||l6|. 
This generates the effective dimension-five neutrino mass operator. The two remaining 
thresholds are treated accordingly, with the number of neutrinos reduced by one each 
time. In practical calculations, this procedure turns out to be computationally demanding. 
We have thus followed the approach of integrating out all the right-handed neutrinos at a 
single intermediate threshold, corresponding to the mass of the lightest of them. We have 
then checked that this approximate treatment does not have any relevant impact on either 
the determination of the GUT-scale parameters or the low-energy predictions of the model. 

Mr, > scale > Mz 

Having constructed the effective theory without right-handed neutrinos at the scale where 
they are integrated out, we use MSSM RGEs [17| to run the gauge couplings, Yukawa 



matrices, soft SUSY breaking parameters and the Wilson coefficient of the neutrino mass 
operator down to the scale Mz- We use two-loop RGEs for dimensionless and one-loop for 
dimensionful parameters. 



So, starting from the fundamental parameters of Step 1, and performing RGE evolution 
through Step 2, one has now the whole set of parameters of the MSSM fixed at Mz or 
below, to tree-level. However, for many quantities, such as gauge couplings, some SUSY 
masses, quark masses, CKM matrix and the Higgs sector, a one-loop determination turns 
out to be mandatory, for different reasons. These issues will be discussed in the next 
section. 

Here we would like to stress a further property of the low-scale MSSM parameters 
implied by the DR model. The off-diagonal entries in the squark mass matrices are gener- 
ated radiatively by the Yukawa couplings; in addition the quark sector features only one 
CP phase^. These two facts allow to classify the low-energy MSSM resulting from the 

^Putting additional phases in the right-handed neutrino mass matrix would communicate corresponding 
phases to the low-energy neutrino sector, without affecting, to a very good approximation, the quark sector. 
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RG evolution of the DR model as belonging to the class of models with minimal flavour 
violation (MFV) |l8|. 

With the above qualifications, using the numerically determined parameters of the 
theory, we can then evaluate effective Hamiltonians for weak decays, in particular for FCNC 
and CP-violating processes. Having these Hamiltonians at hand, it is straightforward to 
evaluate the branching ratios for various low energy processes and to compare them with 
experiments. The latter strategy, allowing to extensively test the model in the flavour 
sector, will be dealt with in Sections ^ and |^ 



4. One-loop improved determination of low-energy observables 

With low-energy observables we denote all the physical quantities that can be directly 
accessed experimentally and compared with predictions of the model. Among these quan- 
tities, some will be used to fix the model parameters, through a fitting procedure, as 
described in Section ^. The others will then be genuine predictions of the model. 

To obtain all physical observables beyond tree- level accuracy at the scale Mz, one needs 
to include one-loop corrections in the relevant formulae defining the observables themselves. 
The inclusion of such corrections is also advocated in the analyses ^. Below, we give 
some details of the procedure in the various cases. 

SUSY masses 

The mass eigenstates for squarks, sleptons, charginos and neutralinos are calculated at the 
weak scale from the tree- level mass matrices |19]. The gluino pole mass is calculated at 
the one-loop level. 

Gauge sector 

We calculate the gauge couplings as{Mz) and acm by including the threshold corrections 
given in |^^. We use tadpole corrections to the Higgs potential |2C, 21 1 to obtain the 
one-loop VEV and include one-loop SUSY corrections to the Fermi constant measured in 
muon decay, G^, as well as to the W and Z boson pole masses [21|. 



Higgs spectrum 

The Higgs spectrum is represented by the masses M^^ , , Mfj+ , Ma of the corresponding 
physical particles. In this model Hq, and A are typically nearly degenerate, with only 
the lightest Higgs mass M^q lying at a substantially lower value. The scale of the masses 
Mhq 1 Mff+ , Ma is set in a non-trivial way by the interplay among fi, tan mn^ ^ and niiQ 
in the equations minimizing the Higgs potential. 

We calculate the pseudoscalar Higgs pole mass Ma following pO] and use it as an 



input to FeynHiggs 2.5.1 |2^, |2^, |2^, 25], which accurately calculates the masses of the 
remaining Higgs mass eigenstates, using Yukawa matrices and soft terms at the EW scale, 
which in turn are the result of the RG analysis. 

We explicitly note that Ma (and with it the other heavy Higgs masses) is typically 
'pushed up' by the upper bound on the BR(i?s ^'^ ^~), so that the Ma value obtained 
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with any given choice of the rest of the parameters can be considered as a lower mass bound 
on the heavy Higgs spectrum. 



Fermion masses and CKM matrix 

At large tan /3, quark masses undergo tan P enhanced corrections that have to be included 
to the tree-level determination, represented by the running Yukawa couplings at Mz |26|. 



These corrections also modify the relations between the original CKM matrix appearing 
in the MSSM Feynman rules and the effective CKM matrix measured in tree-level decays 



1 2^, . We closely follow the line of argument of |£8[ and calculate the one-loop threshold 
corrections to quark and charged lepton mass matrices at Mz, but take into account both 
SUSY and electroweak contributions. After applying the threshold corrections, we use 
three-loop QCD and one-loop QED RGEs to run the five light quark and three charged 
lepton masses down to their respective scales. 

The neutrino masses and the PMNS mixing matrix, on the other hand, are left at the 
tree-level. In fact we find the threshold corrections to these quantities, discussed in p9|| , 
to be numerically negligible. 

5. Basic formulae for FC observables 

In this section we collect formulae for various branching ratios that we will use in our 
numerical analysis. In certain cases we show only the leading contribution for large tan/3, 
in order to provide an intuitive picture of the behaviour. For example, in the case of 
Bs,d /"^/^~ decays, this corresponds to the contribution of Higgs penguins. In the 
actual numerical analysis we include all the relevant contributions, i.e., besides the SM 
one, contributions from charginos, charged Higgses and gluinos. Neutralino contributions 
are generally negligible. 

5.1 Bs^d tJ^~^tJ'~ 

In the SM, the usual Z-penguin and box diagrams result in strongly suppressed branching 
ratios that are sensitive functions of the weak decay constants Fb^ and Fb^ ■ Eliminating 



this dependence with the help of the well- measured mass differences AM^^d [30|, one finds 

m 



BR{B, /U+/i~)sM = (3.37 ± 0.31) x 10"^ , (5.1) 
BR{Bd /U+/x")sM = (1-02 ± 0.09) x 10"^° . (5.2) 



These values should be compared with the present 95% C.L. upper bounds from CDF |l32 | 



BR(S, ^^+/x~)exp < 1.0 X 10-^ BR(5d ^ /x+;u-)exp < 3.0 X 10^\ (5.3) 

that leave still a large room for NP contributions. 

In the MSSM with large tan/3, the helicity suppression in (|5.1| ) and (|5.2| ) is lifted by 
Higgs- mediated neutral currents |33, 34]. Their contributions can be summarized by the 
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approximate formula 25 , |35| 



BR{Bs ^ 3.5 X 10" 



where rrit = mt{fJ.t) and 



tan P 
50 



1.5 ps 



230 MeV 



— 4 



(167r2ey)2 



(1 + es tan /?)2 (1 + eo tan /3)2 



0.040 



(5.4) 



£3 = eo + Vt 



(5.5) 



with eo and ey standing for gluino loop and chargino loop factors, whose full expressions 
can be found in [^ ]. 

With a similar expression for BK{Bd fi'^fi^), one also gets 







2 


\\ytd\] 


2 - 




Fb. 




\yts\ 





Mb, 
Mb. 



n 5 



(5.6) 



where non-leading contributions have been neglected. 

Observing that the ratio ( ^.6P is roughly a factor of ten smaller than the corresponding 
ratio of experimental bounds ( |5.3| ), it is clear that, in our framework, the current BR(i?rf — > 
fi^fj,^) constraint is completely marginal with respect to the Bg counterpart, which is the 
only channel considered in the rest of the analysis. 

We explicitly note that, in eq. (5^), as throughout the text, Vu denotes elements 
of the physical CKM matrix^, to be distinguished for large tan/? from the corresponding 
matrix appearing at the Lagrangian level. The latter is denoted as 'bare' since it does 
not yet include the large tan /?-resummed effects. A similar comment applies to the quark 
masses. In actual numerical calculations, the differences between the physical and the 
'bare' parameters have been taken into proper account as discussed in section ^. 

At this stage, it suffices to state that eo and are at most 0(10^^). Still, with 
tan/3 ~ 50, as characteristic for the DR model, the tan /3-resummed corrections in the last 
factor in (5.4) can be significant and, depending on the sign of fi, can provide an additional 
enhancement of BR(i?s |U^//~) or respectively some suppression relative to the leading 
behaviour (tan/?)^. One finds in the full space of parameters considered 



sign(eo) = sign(e3) = sign(^) 



(5.7) 



and more explicitly. 



2Q;s ^^ tt i 2 I 2 2 / 2\ 

- ^ — ^2 (m^y , ) 



1 At 



167r2 



-H' 



1^ 



(5.8) 
(5.9) 



^In the notation of |28], such matrix elements are written as Vti 



-9- 



where m^, (m^,) are the masses of the i^^ sbottom (stop), At is the soft SUSY breaking 
stop trihnear parameter^ and rrig the gluino mass. The function H2 is defined as 
also IsF 



(see 



H2ix,y) 



X log X 



+ 



ylogy 



{l-x){x-y) {l-y){y-x) 



(5.10) 



We emphasize that in the DR model the parameters entering eo and ey are strongly 
correlated with each other and consequently the range of values which eo and ey can take 
is significantly smaller than in the usual studies of the above formulae, that can be found 
in the literature. 

We also emphasize that in the numerical analysis one can replace the above branching 
ratio with the quantity pCl] 



BR(^, ^ /i+^-) 
AM, 



(5.11) 



thereby eliminating Fb^ that is still inaccurately known. The only hadronic uncertainties 
in the ratio in ( 5.11| ) are present in the non-perturbative factors Bl, that enter AMg only 



linearly and are better known from lattice calculations than the decay constant |37| 



5.2 AM,,d 

The mass differences in the Bg^d — Bg^d systems, AMg 
tan /3, of the following contributions 



consist, in the MSSM at large 



AM,,, = AM|^ + AM,^; + AM,*; + AM|, + AMfj" + AM^^; + AM^^ , (5.12) 

i.e. box diagrams with the SM contribution, with charged Higgses, charginos, gluinos, 
gluino-neutralino and neutralinos and finally neutral Higgs double-penguins, respectively. 
Explicit formulae for AM,,, that include all the important contributions are given in p8[.^ 
The values of i?* are taken from |37]. Besides the dominant SM contribution, the most 

important NP contributions in the DR model are AM^^ and especially AM^J. The latter 



is strictly negative 

AM^P 

mh{^.t) 
3.0 GeV 



= -12.0ps-^ 
rrisifit) 



tan P 
50 



230 MeV 



Wt 



ts\ 



0.040 



0.06 GeV 



(levr^ey) 



(1 + e3tan/3)2(l + eotan/3)2 



(5.13) 



*Our sign convention for At is such that the off-diagonal entry of the tree-level stop mass matrix reads 
mt{At — jj.cot/3). This agrees with the sign convention for Ao in [0, ^ but disagrees with the convention 
used in g 

^We mention that gluino and neutralino box contributions were not considered in [Q. In fact, we find 
these contributions to be very small, but still included them in our numerical analysis. We did not include, 
instead, subleading effects in the Higgs propagator appearing in DP diagrams. The latter have been recently 
addressed in [BSl and play an insignificant role in the present analysis. 
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We recall that, experimentally [|39| , 

(AM,)exp = (17.77 ± 0.10 ± 0.07)/ps , 

to be compared with the UTfit and CKMfitter SM predictions g^, |l[ 

(AM,)SMgt = (18.6 ± 2.3)/ps , (AM,)SM^fitter = (18.9l^|)/ps . (5.14) 

The CKMfitter result still allows for very sizeable NP contributions, while the UTfit 
result bounds |AMs|^^ to be below ~ 3/ps provided other NP contributions in ( ^.12 ) can be 



neglected. We will see in Section |7| that in the DR model AM^ is slightly suppressed relative 
to the SM, but this suppression amounts to at most 5%, in accordance with experimental 
findings. 

Concerning AM^^ , it is suppressed by at least two orders of magnitude relative to 
AM^^ due to rrid/ms and |VjdP/|Vfsp factors and consequently in the DR model AM^ is 
SM-like to a very good accuracy. 



Finally, we note the role played in the present model by the strong correlation |42| 
between the Higgs penguin contributions to BR(i?s^^ and AM^: the enhancement 

of BR(i?s ii^fj,^) in the DR model is correlated with a suppression of AMg. However, 
the data on AMg do not allow for very large enhancements of JiJi{Bs ^ fi^fi^) so that 
both observables turn into constraints on the pseudoscalar Higgs mass Ma- We will come 
back to this point in Section ^. 

5.3 B Xsi 

An important constraint on any NP model is the inclusive decay B — > Xs'y for which the 



data read 1 42 , 



BR{B X^7)cxp = (3.55 ± 0.24 ± 0.10 ± 0.03) x 10"^ , (5.15) 



to be compared with the SM value at the NNLO level I^B] 

BR{B Xsj)sM = (3.15 ± 0.23) x 10^^ . (5.16) 
The inclusion of certain non-perturbative effects decreases this value to (2.98 it 0.26) x 



10 [47|. The latter value, if confirmed, would put some tension on the SM prediction. 



In any case, unless the central experimental value in (5.15) will be significantly decreased. 



NP scenarios predicting BR(i? — > Xgj) to be smaller than the SM value are disfavoured. 
Instead of presenting detailed formulae for B — > in the DR model, which can be 



found in the literature [48, 49, we collect here a number of qualitative properties 

of these formulae that will turn out to be useful in understanding our numerical results. 
These properties are as follows: 

• The charged and neutral Higgs contributions to BR(i? Xs'j) are strictly positive. 

• The sign of the chargino contributions relative to the SM is ruled by the following 
relation 

Cf oc +fiAt tan /3 x sign(C7|^^) , (5.17) 
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with a positive proportionality factor, so it is opposite to that of the SM one for 
// > and At <0 (of. footnote |). 



This shows that the large tan/3 effects in — s- are not as strong as in Bg iJ,~^fj,~, 
where the amplitude behaves as tan^/3/M^. However, they are typically more important 
than in AMs, since in the latter case contributions of Higgs penguins, while behaving as 
tan^/3/M^, are suppressed by the ratio of the external quark masses over M^r. 

Among the NP contributions to B ^ Xg'j, those from charginos are generically the 
largest. In fact, the lightest chargino mass is roughly set by the lowest between /j, and M2 
and in the DR model it turns out to be generically below ~ 200 GeV. On the other hand, 
Higgs contributions are generically small in the DR model, since Mj^ ff+ are pushed up by 
the Bg 1^^ constraint. Consequently, for positive the sign of the Wilson coefficient 
C^ can be reversed relative to Cy^, while AM^, as stated above, cannot be modified by 
more than ~ 5% if the constraint from Bg 1^-^ 1^-^ is taken into account. 

Indeed, for positive fi, in order to pass the B Xgj constraint, the model favours the 
solution Cj{^b) = —Cj^{nb)- We stress here that such solution is a highly conspired one, 
since such equality should hold at the scale, i.e. after running of the coefficients from 
the matching scale. In addition, in this case SUSY is not quite a correction to the SM 
result, but rather the opposite. As a consequence, to address this case, one would need a 
theoretical control on the SUSY part at least as good as that on the pure SM calculation. 
This task is in turn very hard to achieve, since e.g., one would have to accurately know 
where to integrate out the various sectors of the MSSM entering the SUSY contributions 
to — > Xg^. In absence of such knowledge, one can take the approach of matching the 
whole SUSY spectrum at a common reasonable scale. This approach works if SUSY is a 
correction to the SM. But in the present case, the solution Ct{ii},) = —Cj^{^jlij) is extremely 
sensitive to variation on the matching scale, and the theoretical error associated completely 
out of control. We will come back to this point in section ^. 

5.4 B Xsi+£- 

An important observable in our analysis will be the branching ratio for the inclusive decay 
B — > Xsi~^£~ and the related forward-backward asymmetry j4fb- A significant progress 
in calculating this decay and its exclusive counterpart B — >■ K*i~^i~ has been achieved in 
recent years through the calculation of the NNLO QCD corrections. The corresponding 
formulae are very complicated and will not be presented here. They can be found in 

m, n, m, m, n, m, n, n, m. For our discussion it will be sufficient to recall the NLO 



formulae |60, 31 1, keeping only the dipole operator and the operators 

Qq = {sb)v~A{'PfJ')v , Qio = {sb)v^A{'PfJ')A ■ (5.18) 

The contributions of semi-leptonic scalar operators to this decay are much less important 
than in B^^g — > since B Xgi'^i~ is not helicity suppressed. We have taken 



them into account following |62|, but their inclusion in the present discussion would only 



complicate matters without changing the basic picture. 
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Introducing the normalized dilepton mass parameter 



m 



2 ' 



the invariant dilepton mass spectrum in the inclusive decay B 
at NLO as follows 



ds 



sf\VtsfU{s) 



(5.19) 

Xgi'^i' is roughly given 
(5.20) 



where 

U{s) = {l + 2s)(\Cf{s)\' 



\C- 



10 1 



+ 4 1 + 



icf ^ff[2 + uc^^^ MCfis)) , (5.21) 



'7 



with the MSSM expression for the Wilson coefficients C7 and the SM ones for Cg, Cio 



61 



given m [ 

Of particular interest is the (normalized) forward-backward asymmetry in i? — > Xsl^i~ . 
It becomes non-zero only at the NLO level. It is given in this approximation as follows |54| 



ApBis) = -3 Re 



sCf{s) + 2C}; 



C 



10" 



(O)eff 



U{s) 



(5.22) 



The expression for the corresponding asymmetry in the exclusive decay B — > K*i^i~ 
can be found in |58, 59 1. Both asymmetries vanish in the SM at a certain s = sq |66], which 
in the case of the inclusive decay is determined through 



soReCf (so) + 2Cf)^^ 



. 



(5.23) 



In the SM at NLO one finds sq « 0.14. At NNLO the corresponding value is sq = 
0.162 ± 0.008 H, H, 1^, H, g. 

Now, as can be seen in (1^), the very low-s region (s < 1 GeV^) is dominated by 
the coefficient Cj and does not provide more information than already contained in the 
B — > decay. Much more useful is then the low-s region (1 GeV'^ < s < 6GeV'^) which 
is theoretically cleaner than the high-s region, is dominated by the Wilson coefficients Cg 
and Cio and is also very sensitive to the C7-C9 interference in (5.21). For this low-s range 
the world average coming from Belle [39| and BaBar [70| reads 

BR{B )exp = (1.60 ± 0.51) X 10"^ . (5.24) 

Concerning the forward-backward asymmetry ^fBj the only existing positive data 
come from Belle [71|, not yet precise enough to be conclusive on the presence of the zero. 

For our forthcoming discussion, it will be useful to collect the following general prop- 
erties of i? ^ Xgi^i^ in the DR model that will be refined later on: 



'SUSY contributions to C9 and Cio ^ are completely negligible in the DR model. 
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• The Wilson coefficients Cg and Cio receive only small NP contributions so that the 
departures of i? — > Xsi^i~ , ^fb and sq from their SM values are governed by the 
modifications of Cy. 

• Recalling that Cg^ and Cf^ have opposite sign, we observe that the flip of the 
sign of C7 by NP contributions will strongly enhance BR(i? Xsi^i~), while the 
enhancement of [Cyl without the flip of its sign will suppress this branching ratio 
relative to the SM value |5^, |72| 



BR{B Xsi'^r)sM = (1-58 ± 0.10) x lO^*^ 
that is in perfect agreement with experiment. 



(5.25) 



The value of sq is correlated with BR(i? Xg'j) if C7 has the SM sign. It in- 
creases with increasing BR(i? — > ^^7). This is a direct consequence of small NP 



contributions to Cg in most NP models as pointed out in [l73|. 
There is no zero in ApB for sign(C7) = — sign(C7'^). 



5.5 B+ 



T + I/ 



Finally, we consider the tree-level decay 
simply given as follows, 



T~^i/. In the SM, its branching ratio is 



BR(S^ 



T 1^)SM 



GlruB+M^ 
8^^ 



1 



Mi 



m 



B+ 



(5.26) 



As the decay constant -Fg+ ~ Fb^ has still sizeable uncertainties, we consider instead 
the following ratios [^], 



BR{B- 



Svr 



tb+ (AMd)sM 



4 rjB So{mt) Bb^ V 



1 



M2 



Vub 



M2 1 niBi 



tb+{AM, 



s)SM 



^VBSo{mt)BB,M^e ruB, 



Vtd 
2 N 2 



M. 



m 



Vub 



cb 



(5.27) 
(5.28) 



where we used Fb^ ~ Fb+, ttib^ ^ [Vtsj ~ |V"cf,| and the ratio ^ defined as 



Fb.\ Bb 



F, 



(5.29) 



The uncertainties on the right-hand side of ( |5.27| ) and (|5.2^ ) are comparable. In 
( 5.27 ), the only hadronic uncertainty resides in -B_b^, while in ( |5.28| ) there is an additional 
uncertainty in ^. On the other hand, the ratio \Vub/Vcb'^ can be determined from tree-level 
decays without NP pollution, while \Vtd\ can clearly be affected by NP contributions. 

Using the input parameters of Table |2| and |3|, we find from ( 5.27D 



+ + ^ } (0.87 ±0.11) X 10 
\ ;bM I ^^ 2^ _^ Q X 10 



■4 

,-4 



iKbluTfit 
l^wblincl ' 



(5.30) 
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Parameter Value Ref. 


Parameter Value Ref. 


1 n^iT/ , iTTrr^ci sfifin^") iGTii 

|''ub|UTfit o.uuiyLO) \f±<j\\ 

lO^lKblinci 4.49(33) |7| 
W^\Vtd\ 8.49(28) H 
Bb, 1.28(9) 
^ 1.23(6) |77| 
r]B 0.55 


AM fr)s~^l ^7 77(^2^ I75ll 
AMd[ps-i] 0.507(5) |75] 
mt 161.7(2.0) 
ruB, 5.3661(6) \7t\ 
mB+ 5.27913(31) \7t\ 
r5+[10-i2g] 1.638(11) |7e] 



Table 2: Input parameters for the SM prediction of BR(£?+ t^v). Dimensionful quantities are 
expressed in GeV, unless otherwise specified. 



where the first estimate uses the value for \Vub\ resulting from the UTfit analysis of |40|, 
while l^ublinci is the value resulting from inclusive decays alone. The corresponding values 
from ( 5.2^ ) read 



BR(S^ 



(0.82 ±0.12) X lO"'^ 
(1.24 ±0.24) X lO"'^ 



I Kb I UTfit 
iK^'lincl ' 



(5.31) 



showing that the formulae (5.27) and ( |5.28| ) give similar results. 

We observe that the theoretical branching ratio with |T4felinci is closer to the experi- 
mental average between the Belle 



and BaBar [^] results, which reads [pO | 



BR{B^ 



T U 



oxp 



(1.31 ±0.48) X 10" 



(5.32) 



but the large experimental error precludes any clear cut conclusions at present. Yet, simi- 
larly to the case of the B Xg'j decay, extensions of the SM that predict BR(i?"^ — > r"*"/^) 
to be smaller than the SM value seem to be disfavoured at present. 

In this respect two-Higgs-doublet models of type-II, like the MSSM, where each doublet 
couples separately to up- and down-type quarks, are interesting as the interference between 
W and amplitudes is necessarily destructive [Q. One finds then 83 1 



Rbt 



BR{B- 



DR 



BR(S- 



,SM 



1 



m 



B+ 



tan2/3 



m 



H+ 



1 ± eo tan /3 



T/SM 



(5.33) 



We have explicitly shown the dependence on Vub since, in the DR model, the value of \Vub\ 
turns out to be even smaller than \Vub\ extracted from the UT SM fit. 



2)^, where the data seem to be above the SM 



5.6 {g - 2)^ 

In principle we should also consider {g 
expectations by roughly 3a. In many supersymmetric models one finds, for large tan /3, 
/X > and slepton masses 0(400 GeV), additional contributions to {g — 2)^ that allow to 
fit the data. However, in the DR model, the slepton masses are larger than 1 TeV and 
the NP contribution amounts to at most one a of the SM value. Therefore the DR model 
cannot fit the present data on {g — 2)^ and we will not include this observable in the global 
fit, keeping also in mind that the theoretical status of {g — 2)^ is not yet fully satisfactory 
|8|,||,||]. 
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6. General picture 



Having the formulae for the FC observables at hand, we can discuss first the general pattern 
of these observables within the DR model. A detailed numerical analysis will be presented 
in the next section. 

Step 1 

In the DR model, due to the unification of Yukawa couplings, tan /? is forced to be close to 
50. This fact, as already stressed in section requires Ma to be sufficiently large in order 
for the predicted branching ratio of Bs — > in ( |5.4 ) to be consistent with the upper 

bound in ( ^.31 ). Typically we find Ma > 450 GeV. For such large Higgs masses, one has 
Mff+ ~ Ma ~ and this bound is also approximately valid for Mjj+ and Mhq- 

Step 2 

As already stressed in section |5.3| , the large values of Mh+ and Ma imply that NP contribu- 
tions to BR(i? Xg'j) are dominated by charginos, with the positive charged and neutral 
Higgs contributions (as well as those from gluinos and neutralinos) being subleading. 

With /i > it is possible to fit the data on BR(i? — > Xgj) by making the chargino 
contribution so large that Cf^^^ ~ —2 Cy^ at the f^b scale. There are several problems 
with this choice. 



First, there is the problem already stressed at the end of section 5.3. The possibility 
of having Cf^^^ ~ —2 at the fib scale implies that the SUSY contribution is not quite 
a correction to the SM matching condition, but rather the opposite. With such a large 
correction coming from NP, the usual argument of neglecting NLO QCD corrections to 
NP contributions becomes invalid. One would need a control on the NP side at least as 
good as the one present in the SM contribution. This holds not only for the anomalous 



dimension matrix [B7|, but also for the matching conditions of the SUSY contributions 
1 88, 45, Indeed, with the use of only LO matching conditions for SUSY, we find 

a large sensitivity of the finely-tuned condition C^^^^ ~ —2 Cj^ to the choice of the 
matching scale for the SUSY contributions. Such sensitivity becomes even stronger on 
the resulting BR(i? — > Xgj) and this makes a meaningful inclusion of the BR(i? — > Xgj) 
constraint in the numerical analysis practically impossible. 

Second, there is no zero in the forward-backward asymmetry Ap^, which is however 



still elusive experimentally [71|. 



Third and foremost, it has been shown in [91] that this solution for C7 is actually 



excluded by the experimental data on BR(i? — > Xs£^£~), provided the new physics con- 
tributions to the Wilson coefficients Cg^ and Cio are small. In fact, the maximal ranges of 
such contributions in the MSSM with MFV found in ]65t] are too small to bring the theory 
prediction for BR(i? Xgl'^i') in accordance with the experimental data. These findings 
have also been confirmed in f?^ . 

In summary, we want to stress that although we cannot meaningfully take into account 
the case C^^^^ ~ —2Cj^ in the numerical analysis due to the theoretical uncertainties de- 
scribed above, we can exclude this case and impose sign(C7) = sign(C| as a constraint in 
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the global fitting procedure, because of the model-independent arguments brought forward 
in the previous paragraph, which are unaffected by such uncertainties. 



Step 3 

If fi is chosen to be positive and one attempts to be consistent with the data on B ^ Xsi~^i~ 



in (5.24) by keeping Cy to have the SM sign, as discussed in the previous step, the negative 
contribution of charginos tends to suppress BR{B — > Xg'j) below acceptable values and 
can only be tamed by raising the squark masses until these contributions decouple. 

Step 4 

We next move to ^ < 0. In this case, Cj has the same sign as C^^, but due to constructive 
interference between SM, scalar and chargino contributions, it tends to be too large, unless 
squarks are sufficiently heavy. In this respect, we note that, for every given mig, this case 
is not simply a reflection of the corresponding case with > 0. As we will see in the 
numerical section below, for negative /i the lightest squark masses are generically higher 
(> 2 TeV) than in the corresponding positive n case. 

Step 5 

Let us finally look at BR(i?+ — > t'^i'). We have seen that in order to bring the SM value 
for this branching ratio close to its central experimental value, it was necessary to choose 
the tree- level value for \Vub\ > 4 x 10~^. In the DR model, we do not have this freedom as 
\Vub\ is in principle a prediction. In practice it is an outcome of the global fit to the model 
parameters and such fit definitely prefers a low value for |V^fe|, around 3.2 x 10^^. With 
such value and the negative contribution from charged Higgses, we find typically 

BR(B+ T+z^) < 0.6 X 10"^ . (6.1) 



While this is not yet excluded, in view of the large experimental error in ( |5.32 ), also this 



decay could turn out to be problematic for the DR model if the central experimental value 
will remain above 1.0 x 10~^ and the error will decrease by a factor of two. 

Final remarks 

In summary, we have shown that while it is possible through choice of the parameters to 
obtain the agreement of the DR model with a given single observable discussed above, 
simultaneous agreement for all observables is possible at most with very heavy sfermions. 
In the next section, we will present a detailed numerical analysis of these findings. 

7. Numerical analysis 
7.1 Fitting procedure 

We now turn to describe the numerical strategy adopted to test the DR model. As we 



have seen in Section 3.1, the model is completely specified in terms of 24 parameters. 



listed in Table || and here collectively indicated as 'd. After fixing them, it is possible to 
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reconstruct the whole MSSM at low-energy scales, with a well-controlled theoretical error. 
The procedure, based on RGEs, which one adopts to connect the GUT scale model to low 
energies, has been described in Sections |3]^ and ^ 

Once the low-energy MSSM is specified, the model is testable. To this end, one needs 
a suitable set of observables, whose experimental determinations Oi should be as precise as 
possible and, on the theoretical side, calculable within the MSSM with sufficient accuracy. 
Since the MSSM is the low-energy result of the GUT scale model, the theoretical prediction 
for the observable Oi will be functions fi['d] of the model parameters. In order to compare 
theory predictions with experimental values, one defines a suitable x^-function as 



i=l 



where the uncertainty ctj associated with the i observable is defined as 

f^i = \J {(yf)e^p + (c^,^)theo • (7.2) 

Here ((Jj)exp is the experimental RMS error and (o"j)theo an estimate of the theoretical error 
associated with the fi['ff\ calculation. 



The x^-function (7T) is then minimized upon variation of the model parameters 
To this end, we have adopted the minimization algorithm MIGRAD, which is part of the 
CERNlib library |^^. The minimum value for the x^-function provides then a quantitative 
test of the performance of the model in reproducing the observables entering the fit. We 
mention here that, strictly speaking, such test cannot be attached a statistically rigorous 
meaning, i.e. it is not a 'Pearson's test', since, e.g., the x^-entries are not all independently 
measured observables. Nonetheless, the numerical value of the minimum for the function 
( ]7.1| ), as well as the single pulls in its entries, provide a good quantitative indication of the 
detailed performance of the model for the single observables. 

The observables used in the fit are reported in Tables |3[|5|. Concerning the latter, the 
following comments are in order. 

• The observables in Table |^ were already used - among the others - in the previous 
studies 1^, |2| of the DR model, where the very good performance of the model in 
fitting them was demonstrated. We mention that the experimental determination 
of the observables themselves should not rely on any theoretical assumption which 
would be invalidated in the presence of NP, i.e. one should choose observables whose 
determination is NP-independent. This comment applies in particular to CKM- 
related quantities, among which one keeps only those measured through tree-level 
processes and sin 2j3^Xs ' '^hich gives direct access to sin 2/3 since the DR model has 
only one CP phase in the quark sector. 

• The observables in Table on the other hand, represent the real novelty of our study 
with respect to the previous ones. Such FC processes are not calculated after the 
fitting procedure, but instead introduced directly in the x^-function. The procedure 
to calculate these FC observables has been detailed in Section |5|. 
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Observable Value(cTexp) Ref. 



Mw 
Mz 

1/ CKem 

ois{Mz) 
Mt 

mb{mb) 

ms{2 GeV) 
md{2 GeV) 
mu{2 GeV) 



80.403(29) 
91.1876(21) 
1.16637(1) 

137.036 
0.1176(20) 
170.9(1.8) 

4.20(7) 

1.25(9) 
0.095(25) 

0.005(2) 
0.00225(75) 



m 

|76| 
76| 



m 
m 

m 

m 
m 
\m 
m 



Observable 



Value (cJexp) Ref. 



lO^Me 

iKsl 
lO^lKfel 

102|F,fe| 

sin 2/3 
lO^Am^^ 
lO^Amii 
sin^ 2(9i2 
sin^ 2023 



[eV2] 
[eV2] 



1.777(0) 
0.10566(0) 

0.511(0) 
0.2258(14) 

4.1(0.4) 

4.16(7) 
0.675(26) 

2.6(0.2) 
7.90(0.28) 
0.852(32) 
0.996(18) 



m 
m 
m 

HI 

m 
m 
m 
m 
m 



Table 3: Flavour conserving observablcs used in the fit. Dimcnsionful quantities are expressed in 
GeV, unless otherwise specified. 



Observable 



Value(crexp)(o-th e 



Ref. 



lO^ex 

10^ BR(5 ^ Xs-i) 
10^ BR(5 ^ Xsl-^t 
10^ BR(5+ ^ T+v) 
BR(5, ^ /i+/x-) 



€ [1,6] GeV2 



2.229(10)(252) 
35.0(0.4)(3.6) 
3.55(26)(46) 
1.60(51)(40) 
1.31(48)(9) 
< 1.0 X 10-^ 



m 

40| 
111 

m 
s 



Table 4: FC observablcs used in the fit. 



In addition, we included in the fitting function a number of constraints, i.e. those on 
the lightest Higgs mass and on the lightest components of the SUSY spectrum. Table 
^, and the constraint on the BR(i?s fi'^fi^), Table ^. These constraints are in the 
form of suitably smoothened step functions, which are added to the x^-function of 
eq. (|7. l|) . If any of the constraints is violated, the step functions add a large positive 
number to the x^i while for respected constraints the returned value is zero, so that 
the is set back to its 'unbiased' definition (0). 



Further comments on the determination of 
the theoretical errors are in order. First, one 
can note that among the observablcs in Table 
|3[ some have a negligible experimental error. In 
this case, we took as overall uncertainly 0.5% of 
the experimental value, which we consider a re- 
alistic estimate of the numerical error associated 
with the calculations''. Concerning the theoret- 
ical errors on the flavour observablcs (Table H), we note the following: the error on ex is 



Observable 


Lower Bound 


Ref. 


Mho 


114.4 GeV 


m 


mi 


60 GeV 


m 




104 GeV 


m 


m-g 


195 GeV 


m 



Table 5: Mass bounds used in the fit. 



''Note that this error is more conservative than the 0.1% used in ^. 
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Scale 



Mg 



M, 



Mz 



mi, 



MSSMRN RGEs 



Integrate out ur 



MSSM RGEs 



{ Yukawas, SSB, gj, 



SUSY spectrum, 
EWSB conditions, 
one-loop corrections to 
masses and mixings 



^CNC observables [■ 



Fit of 

textures (e, e', . . .) 

Mg, OiG, £3 




Figure 1: Schematic chart of the strategy followed in the fitting procedure. 



basically that on the lattice parameter Bk] the error on AM<j/AM(i keeps into account that 
on the SM contribution, dominated by and that on the NP contributions, dominated by 
the scalar ® Pr operators; the error on BR(i3+ — > t'^v), after normalization as in eq. 
(5.27), is only that on B^; the error on BR(i? ^sl) is taken as twice the total theoreti- 



cal error associated with the SM calculation 1^^; finally the error on BR(i? — > Xsl'^i~) is 
taken as 25% of the experimental result, and is estimated from the spread of the theoretical 
predictions after variations of the scale of matching of the SUSY contributions. 

We next turn to the generic strategy adopted to minimize the x^-function with respect 
to the parameters We note that, among them, ^, and tan (3 are those responsible 
for EW symmetry breaking, and the x^-function manifests a particularly sensitive depen- 
dence on them, especially on mn^^- As a consequence, such parameters are varied first 
(keeping the other fixed to initial guesses), in order to successfully find an EW symme- 
try breaking minimum, thereafter varying the rest of the parameters. This procedure is 
schematically described in the flow-chart of Fig. ||. As a final step, all the parameters are 
varied simultaneously. 

Among the model parameters, the ones in the SUSY sector are particularly interesting 
since they set the scale of the SUSY particles' masses. This is especially true for mig and 
/X. All the other parameters can be left free in the fit, since their typically allowed range of 
variation is quite narrow. In particular, when mig is fixed and jj, is positive, the fit prefers 
regions of the remaining parameter space such that 

IJ,,Mi/2 < mi6 , - ^0 - 2mi6 , (7.3) 
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which is favoured by third generation Yukawa unification |94, 96 1. Concerning the 
first of relations (7.3), is bounded from above because otherwise the bottom mass is 



pushed up beyond acceptable values by large gluino corrections. In fact, we find that M1/2 
is most of the times chosen in the range [140, 400] GeV, where the lower bound results from 
the chargino mass bound in Table |5[ 

The second of relations ( |7.3| ) leads to an inverted scalar mass hierarchy i.e. heavy 
first and second generation sfermions, but lighter third generation sfermions. For the values 
of mi6 considered here, namely mie > 4 TeV, and for ^ > 0, this condition also helps to 
obtain the correct prediction for rrib p^j . 

On the other hand, the allowed interval for fj, is generically wider, for every fixed value 
of miQ. As a consequence, our main strategy is to study the model behavior for different 
choices of {mi6,/i}, and for each of them, let the rest of the parameter space free to be 
determined by the minimization procedure. In the next section we now turn to describe 
the various scenarios considered in the {mi6,/u} plane. 

7.2 Scenarios 

We considered increasing values of mig starting from 4 TeV, which represents the 'mini- 
mum' value for successful fits to the observables of Table |^ |]|, |2|. For each fixed value of 
miQ we then studied the /i dependence by performing fits with different initial guesses for 
this parameter. All input values for the model parameters corresponding to these scenarios 
are listed in Table |6|. In the following, we describe in detail our findings. 

7.2.1 mi6 = 4 TeV, /x > 

Given the inverted scalar mass hierarchy, the relatively low value of m-ig leads to stop 
masses below 1 TeV, resulting in a large chargino contribution to BR(i? — > Xgj)- 

As a matter of fact, the preferred region of 
NP contributions to the latter decay mode re- 
verses the sign for Ci{iJL}y): C7{fib) — ~Cj^{fJ-b)- 
On the fine-tuned character of this case we have 
already commented in Sections 5^ and |6|. We 
stress again that within the DR model, this solu- 
tion is not viable in view of the implied enhance- 
ment to the branching ratio of i? ^ Xgl^i^- 




-0.5 




SM H 



5 SM+NP 



Model- independent analyses |91, 72] show in fact 
that to compensate for such enhancement, one 
would need substantial contributions to the Wil- 
son coefficients Cg^ and Cio of the operators in 
eq. ( [5.18| ). However, within the DR model, these 
Wilson coefficients are always SM-like to an ex- 
cellent approximation. 

As a consequence of the above, we have studied the viability of having the SM sign in 
CjifJ-b) by imposing this condition as a constraint on the function. The typical fit result 
in this case is illustrated in Table |^. 



Figure 2: Contributions to C7{fib) for the 



scenario of section 7.2.1 



- 21 - 



The fit displays the main problem of the model in this regime for {mig,//}, i.e. a 
5(7 discrepancy in the predicted BR(i? ^sl)- Even imposing the SM sign on Cjdib), 
the contribution from charginos is still too large in magnitude, and none of the other NP 
contributions is able to compensate for it. The situation is illustrated in Fig. |2[ 

A second, though less severe, problem is the predicted value for BR(i?+ t'^'^), which 
is roughly 2a too low with respect to the experimental average ( |5.32| ). This problem is 
strictly connected to the quite low value for \Vub\ ~ 3.2 x 10~^ predicted by the model. We 
found this feature to hold irrespective of the values chosen for {mi6,/x}, so that it should 
be connected to the specific Yukawa textures of the model. 

7.2.2 mi6 = 6 TeV, /x > 

Since the BR(i? Xg^) problem is related to the low value for fi required (at least for 
positive /i) by mig = 4 TeV, we have tried to increase the latter in order to understand 
how fastly decoupling is effective in mildening the problem. 

For mi6 = 6 TeV and n > the preferred 
range for ^ is roughly ^ G [800, 1000] GeV. A 
typical fit result is displayed in Table ^, corre- 
sponding to the case /i = 953 GeV. As a matter 
of fact, the discrepancy in B ^ Xgj is tamed 
to roughly 2.3a, since chargino contributions are 
less important than in the miQ = 4 TeV cases. 
The various contributions to Cid^h) for the fit in 
Table |8| are displayed in Fig. |3|. By decreasing /j, 
below « 800 GeV, the prediction for B Xsj 
fastly worsens. For example, for a converged fit 
with n = 430 GeV, we found that the discrep- 
ancy is already at the 4.2(T level. Concerning 
B^ T~^i>, as anticipated above, the predicted 
rate remains always roughly 2a off. 

As a final remark, while the case mig = 6 TeV allows, for suitable /i, a smaller discrep- 
ancy in 5 —J- X<j7, the latter comes with the price of a much higher mass for the lightest 




SM+NP 



Figure 3: Contributions to C-;{^h) for the 
scenario of section 7.2.2| . 



up- type squark, as evident by comparing the corresponding values in Table 11 



7.2.3 mi6 = 10 TeV, /i > 

Increasing mig further, the chargino contribution to BR(S Xg'^) becomes comparable 
in size to the charged Higgs contribution, resulting in acceptable values for this branching 
ratio. For example, in the fit of Table ^, the pull in this observable is reduced to 1.3(7. 
On the other hand, the pull resulting from B^ — > t~^v is not significantly ameliorated 
compared to the previous cases. In addition, the lightest squark is as heavy as 1.9 TeV. 



7.2.4 mi6 = 4 TeV, /i < 



We explored also the case with negative /i. In this instance, relations in eq. (7^) (with 
fx — > |//|) do not apparently need to be fulfilled. As a matter of fact, we found satisfactory 
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fits for a quite wide range of /i: ^ G —[2100, 400] and Aq is always lower in magnitude than 
the value required by the second relation in eq. ( [7.3[ ), typically leading to very small At. 
As a consequence, the squark mass spectrum does not fulfill an inverted mass hierarchy 
1 97, 96 1 and squark masses are generically very heavy. In this case, also heavy Higgses are 
found to have generically larger masses, ^1-5 TeV. A typical result is shown in Table 10 
and the displayed features remain basically the same in the full allowed range for /i. 

We observe that, in this case, negative values of /i, large squark masses and small values 
for At imply small threshold corrections to and therefore allow to have successful Yukawa 
unification away from the inverted mass hierarchy condition. We have then investigated 
whether an acceptable fit away from inverted mass hierarchy could also be obtained for 
/i > 0^. In the latter case, the most important contributions to the mf, threshold correction 
have the same sign and consequently one generically needs larger squark masses than for 
H < 0, in order to reproduce the right rrit, value. In fact we find that, unless mig > 6 TeV, 
the prediction on mi, is Aa too large^ and consequently fits with /i > and away from the 
inverted mass hierarchy perform worse than the corresponding negative /i cases. 

We also note that, in the {r?T,i6,/x} mass scenario considered in the present subsection, 
the predicted value for BR(i? Xg'j) is always larger than the SM prediction and close to 
the experimental value. In fact, small At implies negligible chargino contributions, so that 
the main correction is the one from Higgses. However, since the lightest stop is around 
2.6 TeV, this scenario clashes with the motivation for SUSY as a solution to the Higgs 
fine-tuning problem. 
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7.3 Results 

Considering the discussion in the last sec- 
tion, it is apparent that for positive 
the tension between the three decays B 
Xsj, B Xsi^i^ and Bg can be 

relieved by raising the universal sfermion 
mass miQ beyond 8 TeV. This is demon- 
strated in Fig. ^, showing the total 
and the contribution of i? ^ Xg^ to the 
for increasing values of mig. For neg- 
ative fi, fits with comparable can be 
achieved with lower values of mig, but not 
with lighter squarks, as the lightest squark 
is still very heavy in these cases. 

However, it is well-known that the su- 
persymmetric solution to the gauge hier- 
archy problem requires light third genera- 
tion sfermions. Therefore, light stops are 
favourable to reduce fine-tuning. To show 
the amount of splitting between fermion and sfermion masses needed in the DR model, we 
plot the lightest stop mass m^^ versus the total ^or all fits with positive as well as neg- 
*We warmly thank R. Dermi'sek for drawing this point to our attention. 

^One should also take into account the quite precisely known value for rrib assumed in the present paper. 



4 TeV, 6TeV, 8 TeV, lOTeV, 4 TeV, 
378 GeV 953 GeV 1360 GeV 1200 GeV-2106 GeV 

{mi6 . /"I 

Figure 4: Total (whole bar) and contribution 
of S — > Xsj to the (red bar) for the best fits 
with positive /i and mie = 4, 6, 8 and 10 TeV, 
respectively, as well as the best fit with negative 
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ative n we obtained (see Fig. ^). There is obviously a strong correlation between the stop 
mass and the quality of the fit, demonstrating that m^^ has to be at least as large as 1.8 TeV 
to cure the problems with the three aforementioned decays. This is significantly heavier 
than the masses considered in [|l|, ^ and may be difficult to reconcile with naturalness. 

In addition, a number of problems of 
the model persist even for very large mig. 
These are the issues related to the Yukawa 
textures: the up-quark mass, Vub and sin 2/3, 
and as a result also AMs/ AM^ and —>■ 
T'^v. This is why there are no points with 

^ 25 in Fig. ^. The reason for the 
much higher contribution from flavour 
conserving quantities as compared to ||l|, ^ 
is mainly due to updated experimental val- 
ues and reduced experimental errors, espe- 
cially in \Vus\, \Vub\ and m^. 
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Figure 5: Total vs. the lightest stop mass for 
all obtained fits. Red circular points are fits with 
positive jjL, blue squares with negative /j,. 



For a given value of mig, we found 
that successful fits could only be obtained 
in a limited range of allowed values for 
The preferred value of ^ increased with in- 
creasing mi6, as can be seen for the best 
fit points in Fig. ^. This fact could give 

rise to an additional problem. Although we did not include the dark matter density as a 
constraint in our analysis, we conjecture that such large values of n would give rise to a 
relic abundance of neutralinos incompatible with the WMAP measurements, for the follow- 
ing reason. Because of the extremely heavy sfermions, the dominant annihilation channel 
for the neutralino in this scenario is through an s-channel pseudoscalar Higgs. While this 
rate is already suppressed by the large Ma, the coupling of neutralinos to the pseudoscalar 
Higgs is additionally suppressed by large ^u. This would result in an overabundance of relic 
neutralinos incompatible with observations. Solving this problem by resonant neutralino 
annihilation with m^^o ~ M^/2, as was advocated in an extensive analysis of dark matter 
in this class of models |^, |9^, is not possible because of the large Ma and the small -M1/2 
preferred by the fit. 



8. Conclusions 

In this paper we have performed a detailed analysis of the SO(IO) SUSY GUT model with 
family symmetry of Dermisek and Raby ||l|, || . 

This model is entirely specified in terms of 24 parameters. Once they are fixed, the 
whole MSSM parameter space (including its SM subset) is predicted at low energies with 
the help of RGEs. The common dependence on the model parameters strongly correlates 
all the low energy observables, in contrast with the direct consideration of the MSSM at 
the EW scale, where the CKM parameters and the fermion masses are fully independent of 
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the SUSY particle speetrum. We find that the DR model gives a satisfactory description 
of the quark and lepton matrices as well as of the PMNS and CKM mixing matrices 
with one exception: the CKM element \Vub\ turns out to be significantly smaller than 
l^blinci extracted from inclusive tree- level decays and even smaller than |V^b|exci- The 
above findings are mostly a confirmation of previous studies of the model. 

The main novelty of our study, with respect to similar analyses of SUSY GUT models 
found in the literature, is that we analyze simultaneously the mass spectra of quark and 
leptons, the CKM and PMNS mixing matrices, the SUSY mass spectrum and its implied 
corrections to the FC processes Bg A*^A*~, B Xg^, B Xs£~^i~, 5+ t'^u and the 
B^ g — Bd^s mass differences AM^^g- The performance of the model is assessed by means of 
a global fit to the above mentioned observables. 

The inclusion of the FC processes listed above turns out to be a crucial test of the 
mass hierarchies predicted by the model for the SUSY spectrum. In fact, such hierarchies 
unavoidably manifest themselves in loop corrections, and FC observables remain the most 
sensitive probes of such corrections. Our analysis demonstrates that the simultaneous 
description of all the FC processes listed above is a serious challenge for the DR model. 
In view of the specific way this failure is realized, we suspect that this is a problem of a 
wider class of SUSY GUTs in which the presence of Yukawa unification implies tan (3 ~ 50, 
unless non-minimal sources of flavour violation are introduced. 

Our main message is the following one. To really assess the viability of models for 
flavour parameters, it is essential not only to verify their ability to reproduce quark and 
lepton mass spectra and mixing matrices - in itself an already notable achievement - 
but also to test the consistency with the data on available FC processes, since the latter 
have a simultaneous sensitivity to mixing matrices and new particles' spectra. In the DR 
model example, FC processes arc in fact the best probes available to the SUSY part of 
the spectrum, where information from direct detection is missing. It turns out that the 
DR model - otherwise successful for quark and lepton masses as well as for the CKM and 
PMNS matrices - is challenged only when specifically tested in the simultaneous description 
of quark FCNC processes. The failure in the description of these data makes the viability 
of the DR model questionable from the present perspective, but hopefully offers insights 
on further lines of development along similar classes of models. 
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Note added 

During the completion of the present work, a new bound on the branching ratio for Bg 
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/i^/i has been presented at the HEP 2007 conference [|100|| . The latter results from a 
combined analysis of the CDF and D0 data and reads 



BR{Bs /U+^-)cDF+D0 < 5.8 X 10~^ (95% C.L.) . (8.1) 

This bound represents a considerable improvement over the one given in eq. ( ^.3[ ). 

In addition, at the SUSY 2007 conference it was presented a new (preliminary) result 
on the BR(i?~'~ t^u) from the BaBar collaboration |101[| , which reads 

BR(5+ ^ r+l/)BaBar,prelim. = (1-2 ± 0.4stat ± 0.3bkg ± 0.2eff) X 10"^ . (8.2) 

Performing the weighted average between the result in eq. ( |8.2D and the Belle result ||78| 
one obtains 



BR(5+ ^ T+l/)new = (1-41 ± 0.43) X 10" 



(8.3) 



The new averages in eqs. (8.1) and (B^) further strengthen our conclusions. 
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mi6 


A nnn 
4UUU 


DUUU 


1 nnnn 
iUUUU 


A nnn 
4UUU 






yoo 


1 onn 


— ZiUO 


Ml/2 


147.3 


145.6 


146.7 


229.9 




-7787.4 


-11924.0 


-20070.0 


-630.13 


tan/3 


49.9 


48.8 


48.7 


49.3 


l/ac 


24.7 


24.7 


24.7 


24.6 


Mg/10^*5 


2.45 


3.11 


4.8 


4.96 


£3/% 


-3.81 


-4.12 


-5.37 


-5.6 


(mH„/mi6) 


1.59 


1.55 


1.57 


0.52 


/ / \ 9 

(mH^/mie) 


1.86 


1.79 


1.8 


1.0 




1.059 


1.05 


1.072 


f\ f\C\f\ 

0.8999 




-74.85 


-70.51 


-71.93 


-66.94 




3244.0 


3053.0 


3069.0 


2718.0 


A 


0.618 


0.583 


0.582 


0.578 


e 


0.0473 


0.048 


0.0477 


0.048 


e' 


-0.0034 


-0.00338 


-0.00342 


-0.00356 


IpI 


0.0566 


0.0584 


0.0574 


0.0567 


cr 


1.02 


1.03 


1.03 


0.987 


e 


0.00964 


0.00957 


0.00967 


0.00999 


lel 


0.134 


0.134 


0.134 


0.14 


arg /9 


3.93 


3.92 


3.94 


3.94 


arg a 


0.641 


0.617 


0.622 


0.751 


arg e 


0.484 


0.492 


0.491 


0.505 


arg^ 


3.61 


3.6 


3.61 


3.58 



Table 6: Input parameters for the fits presented in section (cf. Table |l|). Dimensionful 
quantities are given in units of GeV. 



-32- 



Observable 


Exp. value 


Fit value 


Pull {a) 


Mw 


80.403 


80.4 


0.1 


Mz 


91.1876 


90.6 


1.3 


Gf X 10^ 


1.16637 


1.16 


0.7 


l/Oem 


137.036 


136.4 


0.9 


as{Mz) 


0.1176 


0.115 


1.1 


Mt 


170.9 


171.4 


0.2 


nibimb) 


4.2 


4.31 


1.5 


mc{mb) 


1.25 


1.15 


1.2 


m^(2GeV) 


0.095 


0.107 


0.5 


md(2GeV) 


0.005 


0.00741 


1.2 


m„(2GeV) 


0.00225 


0.00462 


3.2 


Mr 


1.777 


1.77 


0.4 




0.10566 


0.106 


0.1 


Me 


0.000511 


0.000511 


0.0 


\Vus\ 


0.2258 


0.225 


0.6 


\Vub\x 


4.1 


3.26 


2.1 


\Vcb\ 


0.0416 


0.0417 


0.1 


sin 2/3 


0.675 


0.637 


1.4 


Amii X lO^i 


2.6 


2.6 


0.0 


Aml^ X 10^3 


7.9 


7.9 


0.0 


sin^ 2012 


0.852 


0.85 


0.1 


sin^ 2^23 


0.996 


1.0 


0.2 


CK X 10^ 


2.229 


2.32 


0.4 


BR{B Xs-f) X 10^ 


3.55 


0.885 


5.0 


BR{B Xsl-^r ) X 10*^ 


1.6 


1.8 


0.3 




35.05 


29.8 


1.4 


BR(5+ ^ r+i/) X 10^ 


1.31 


0.336 


2.0 






total x^: 


58.3 



Table 7: Fit results for the case mie = 4 TeV, n = 378 GeV. The pull for the i"^ observable 
represents the square root of the corresponding entry in the function (^). Corresponding 
predictions are reported in Table O. Dimcnsionful quantities are given in units of GeV. 
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Observable 


Exp. value 


Fit value 


Pull {a) 


Mw 


80.403 


80.5 


0.2 


My 


91.1876 


90.6 


1.2 


Gy X 10^ 


1.16637 


1.16 


0.5 


1 / OLr,m 


137.036 


136.5 


0.8 


\ Zj J 


0.1176 


0.116 


0.5 


Mt 


170.9 


169.8 


0.6 


TYlui TTlh) 

" ^0\" ^0 J 


4.2 


4.29 


1.3 




1.25 


1.14 


1.2 


ms{2GeV) 


0.095 


0.106 


0.4 


md(2GeV) 


0.005 


0.00727 


1.1 


m„(2GeV) 


0.00225 


0.00465 


3.2 


Mr 


1.777 


1.77 


0.3 




0.10566 


0.106 


0.1 


Me 


0.000511 


0.000511 


0.0 


\Vus\ 


0.2258 


0.225 


0.6 


\Vub\x 10^ 


4.1 


3.26 


2.1 


\Vcb\ 


0.0416 


0.0417 


0.1 


sin 2/3 


0.675 


0.638 


1.4 


Aml^ X 10^^ 


2.6 


2.6 


0.0 


Am|i X 10^3 


7.9 


7.9 


0.0 


sin^ 26*12 


0.852 


0.852 


0.0 


sin^ 26*23 


0.996 


0.997 


0.1 


ex X 10^ 


2.229 


2.31 


0.3 


BR{B Xsj) X lO-^ 


3.55 


2.34 


2.3 


BR{B Xsi+i-) X 10*^ 


1.6 


1.62 


0.0 


AM,/AMd 


35.05 


30.0 


1.4 


BR(5+ ^ T+v) X lO'^ 


1.31 


0.398 


1.9 






total x^'- 


35.6 



Table 8: Fit results for the case mie = 6 TeV, /i = 953 GeV. Corresponding predictions are 
reported in Table ^ Dimensionful quantities are given in units of GeV. 
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Observable 


Exp. value 


Fit value 


Pull {a) 


Mw 


80.403 


80.6 


0.5 


My 


91.1876 


90.7 


1.1 


Gy X 10^ 


1.16637 


1.16 


0.3 


1 / OLr,m 


137.036 


136.8 


0.4 


\ Zj J 


0.1176 


0.117 


0.2 


Mt 


170.9 


170.6 


0.2 


TYlui TTlh) 

" ^0\" ^0 J 


4.2 


4.22 


0.3 




1.25 


1.14 


1.2 


ms{2GeV) 


0.095 


0.107 


0.5 


md(2GeV) 


0.005 


0.00741 


1.2 


m„(2GeV) 


0.00225 


0.00461 


3.1 


Mr 


1.777 


1.78 


0.1 




0.10566 


0.106 


0.1 


Me 


0.000511 


0.000511 


0.0 


\Vus\ 


0.2258 


0.225 


0.6 


\Vub\x 10^ 


4.1 


3.26 


2.1 


\Vcb\ 


0.0416 


0.0416 


0.1 


sin 2/3 


0.675 


0.639 


1.4 


Aml^ X 10^^ 


2.6 


2.6 


0.0 


Am|i X 10^3 


7.9 


7.9 


0.0 


sin^ 26*12 


0.852 


0.852 


0.0 


sin^ 26*23 


0.996 


1.0 


0.2 


ex X 10^ 


2.229 


2.33 


0.4 


BR{B Xsj) X lO-^ 


3.55 


2.86 


1.3 


BR{B Xsi+i-) X 10*^ 


1.6 


1.62 


0.0 


AM,/AMd 


35.05 


31.1 


1.1 


BR(5+ ^ T+v) X lO'^ 


1.31 


0.517 


1.7 






total x^'- 


27.4 



Table 9: Fit results for the case mig — 10 TeV, /i = 1200 GeV. Corresponding predictions are 
reported in Tabic ^ Dimcnsionful quantities are given in units of GeV. 



-35- 



Observable 


Exp. value 


Fit value 


Pull {a) 


Mw 


80.403 


80.7 


0.7 


My 


91.1876 


90.7 


1.1 


Gy X 10^ 


1.16637 


1.17 


0.2 


1 / OLr,m 


137.036 


136.9 


0.3 


\ Zj J 


0.1176 


0.118 


0.0 


Mt 


170.9 


170.5 


0.2 


TYlui TTlh) 

" ^0\" ^0 J 


4.2 


4.19 


0.1 




1.25 


1.14 


1.2 


ms{2GeV) 


0.095 


0.0999 


0.2 


mw(2GeV') 


0.005 


0.00716 


1.1 


m„(2GeV) 


0.00225 


0.00446 


3.0 


Mr 


1.777 


1.78 


0.1 




0.10566 


0.106 


0.2 


Me 


0.000511 


0.000511 


0.1 


\Vus\ 


0.2258 


0.224 


1.2 


\Vub\x 10^ 


4.1 


3.26 


2.1 


\Vcb\ 


0.0416 


0.0416 


0.0 


sin 2/3 


0.675 


0.64 


1.3 


Aml^ X 10^^ 


2.6 


2.6 


0.0 


Am|i X 10^3 


7.9 


7.9 


0.0 


sin^ 26*12 


0.852 


0.851 


0.0 


sin^ 26*23 


0.996 


0.996 


0.0 


ex X 10^ 


2.229 


2.35 


0.5 


BR{B Xsj) X lO-^ 


3.55 


3.34 


0.4 


BR{B Xsi+i-) X 10*^ 


1.6 


1.63 


0.0 


AM,/AMd 


35.05 


31.4 


1.0 


BR(5+ ^ r+i/) X lO'^ 


1.31 


0.59 


1.5 






total x^'- 


24.5 



Table 10: Fit results for the case mie = 4 TeV, ^ = —2106 GeV. Corresponding predictions are 
reported in Table ^ Dimcnsionful quantities are given in units of GeV. 
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4000 


6000 


10000 


4000 




378 


953 


1200 


-2106 


BR{Bs /ti+/U-) X 10^ 


8.6 


7.7 


2.1 


0.33 


So 


0.022 


0.13 


0.14 


0.16 


BR(^i ^ 67) X 10^3 


0.36 


0.021 


0.0026 


0.011 




+5.8 


+1.6 


+0.52 


-2.9 




126 


129 


129 


119 


Ma 


507 


559 


842 


1800 




640 


1172 


1903 


2627 




895 


1475 


2366 


2488 




1510 


2419 


3933 


2931 




60 


60 


60 


94 




115 


119 


120 


189 




462 


478 


506 


703 



Table 11: Predictions for the scenarios presented in section 7.2. Masses are given in units of GeV. 
So is the zero position of the forward-backward asymmetry in _B ^ Xgl^t^ (cf. Section |3.4[ ) . The 
quantity 5a^^^^ is the SUSY contribution to = [g — 2)^/2, which is currently measured to be 
+27.6 X 10"^" larger than the SM prediction |85, 86 , with an uncertainty of about 8 x 10' 
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